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ABSTRACT 

We are presenting a detailed parameter study of the time-dependent electron 
injection and kinematics and the self-consistent radiation transport in jets of in- 
termediate and low-frequency peaked BL Lac objects. Using a time-dependent, 
combined synchrotron-self-Compton and external-Compton jet model, we study 
the influence of variations of several essential model parameters, such as the elec- 
tron injection compactness, the relative contribution of synchrotron to external 
soft photons to the soft photon compactness, the electron-injection spectral in- 
dex, and the details of the time profiles of the electron injection episodes giving 
rise to flaring activity. In the analysis of our results, we focus on the expected 
X-ray spectral variability signatures in a region of parameter space particularly 
well suited to reproduce the broadband spectral energy distributions of inter- 
mediate and low-frequency peaked BL Lac objects. We demonstrate that SSC- 
and external-Compton dominated models for the 7 -ray emission from blazars are 
producing significantly different signatures in the X-ray variability, in particular 
in the soft X-ray light curves and the spectral hysteresis at soft X-ray energies, 
which can be used as a powerful diagnostic to unveil the nature of the high-energy 
emission from BL Lac objects. 
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Subject headings: galaxies: active — galaxies: jets — radiation mechanisms: 
nonthermal 


1. Introduction 

The class of objects referred to as blazars consists of the most extreme examples of 
active galactic nuclei (AGNs), namely 7-ray loud, flat-spectrum radio quasars (FSRQs), and 
BL Lac objects. They have been observed in all wavelength bands — from radio through 
very-high energy (VHE) 7-ray frequencies. More than 65 blazars have been identified as 
sources of > 100 MeV emission detected by the EGRET telescope on board the Compton 
Gamma-Ray Observatory (CGRO) (Hartman et al. 1999), and at least 5 blazars have now 
been detected at VHE 7-rays (> 350 GeV) by ground-based air Cerenkov telescopes (for 
a recent review see, e.g., Buckley (2001)). Blazars exhibit variability at all wavelengths on 
time scales — in some cases — down to less than an hour (Gaidos et al. 1996). 

The broadband continuum spectra of blazars are dominated by non-thermal emission 
and consist of at least two clearly distinct, broad spectral components. A sequence of sub- 
classes of blazars can be defined through the peak frequencies and relative uF u peak fluxes of 
those components, which also appear to be correlated with the overall bolometric luminosity 
of the sources (Fossati et al. 1998): In the case of flat-spectrum radio quasars (FSRQs), 
the low-frequency (synchrotron) component extends from radio to optical/UV frequencies, 
with a peak frequency generally in the mm or IR band; the high-frequency component 
extends from X-rays through GeV 7-ray energies, with a vF„ peak frequency corresponding 
to ~ 10 MeV 1 GeV. No FSRQ has so far been detected by ground-based air Cerenkov 
telescope facilities at energies > 100 GeV, although in flaring states the 7-ray vF v peak flux 
of FSRQs dominates over the low-frequency emission by up to ~ 1 order of magnitude. In the 
case of high-frequency peaked BL Lac objects (HBLs), the low-frequency component often 
extends far into the X-rays, with peak frequencies ranging from the UV/soft X-ray to the hard 
X-ray regime, depending on the source and its state of activity; the high-energy component 
of HBLs extends from hard X-rays into the VHE 7-ray regime. All blazars detected at VHE 
7-ray energies to date are HBLs. In spite of extending to extremely high photon energies, 
the vF v peak flux of the 7-ray component of HBLs is generally at most comparable to 
the spectral output in the low-frequency component. In terms of their overall bolometric 
luminosity, FSRQs appear to be several orders of magnitude more powerful sources than 


Chandra Fellow 



- 3 - 

HBLs. Apparently intermediate between the FSRQs and the HBLs are the low-frequency 
peaked BL Lac objects (LBLs). The peak of their low-frequency component is typically 
located at IR or optical wavelengths, their high-frequency component peaks typically at 
~ several GeV, and the 7 -ray output is of the order of or slightly higher than the level of 
the low-frequency emission. 

Although all extragalactic sources detected by ground-based air Cerenkov telescope fa- 
cilities to date are HBLs, the steadily improving flux sensitivities and decreasing energy 
thresholds of those instruments provide a growing potential to extend their blazar source 
list towards intermediate and even low-frequency peaked BL Lac objects. The detection of 
such objects at energies ~ 40 — 100 GeV might provide an opportunity to probe the intrinsic 
high-energy cutoff of their spectral energy distributions (SEDs) since at those energies, 77 
absorption due to the intergalactic infrared background is expected to be negligible at red- 
shifts of z < 0.2 (de Jager k Stecker 2002). Such detections should significantly further our 
understanding of the relevant radiation mechanism responsible for the high-energy emission 
of blazars and the underlying particle acceleration mechanisms. 

The low-energy component of blazar SEDs is well understood as synchrotron emission 
from ultrarelativistic electrons in a relativistic jet directed at a small angle with respect to 
the line of sight. Several scenarios have been proposed concerning the acceleration of such 
ultrarelativistic electrons, including impulsive injection near the base of the jet (Dermer k 
Schlickeiser 1993; Bottcher, Mause, k Schlickeiser 1997), individual shock waves propagating 
along the jet (Marscher k Gear 1985), or internal shocks from the collisions of multiple shells 
of material ejected into the jet structure (Spada et al. 2001). In the framework of leptonic 
models, high-energy emission will result from Compton scattering of lower- frequency photons 
off the relativistic electrons. Possible target photon fields for Compton scattering are the 
synchrotron photons produced within the jet (the SSC process; Marscher k Gear (1985), 
Maraschi, Ghisellini, k Celotti (1992); Bloom k Marscher (1996)), or external photons (the 
EC process). Sources of external seed photons include the UV - soft X-ray emission from the 
disk — either entering the jet directly (Dermer, Schlickeiser, k Mastichiadis 1992; Dermer k 
Schlickeiser 1993) or after reprocessing in the broad line region (BLR) or other circumnuclear 
material (Sikora, Begelman, k Rees 1994; Blandford k Levinson 1995; Dermer, Sturner, k 
Schlickeiser 1997) — , jet synchrotron radiation reflected at the BLR (Ghisellini k Madau 
1996; Bednarek 1998; Bottcher k Dermer 1998), or the infrared emission from circumnuclear 
dust (Blazejowski et al. 2000; Arbeiter, Pohl, k Schlickeiser 2002). 

It has been suggested that the sequence of spectral properties of blazars from HBLs via 
LBLs to FSRQs can be interpreted in terms of an increasing total power input into non- 
thermal electrons in the jet, accompanied by an increasing contribution of external photons 
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to the seed photon field for Compton upscattering (Madejski 1998; Ghisellini et al. 1998). 
Detailed modeling of blazars in the different sub-classes (FSRQs, LBLs and HBLs) seems 
to confirm this conjecture: HBLs are generally well fitted with pure SSC models, while 
FSRQs often require a dominant contribution from external photons to the soft seed photon 
field in order to produce the high-energy emission (for a recent review, see, e.g., Bottcher 
(2002)). A moderate contribution of external soft photons appears to be required in LBLs 
like BL Lacertae (Madejski et al. 1999; Bottcher & Bloom 2000) or W Comae (Bottcher, 
Mukherjee, & Reimer 2002). 

As mentioned earlier, blazars tend to exhibit rapid flux and spectral variability. The 
variability is most dramatic and occurs on the shortest time scales at the high-energy ends 
of the two nonthermal spectral components of their broadband SEDs. Particularly interest- 
ing variability patterns could be observed at X-ray energies for those blazars whose X-ray 
emission is dominated by synchrotron emission. Observational studies of X-ray variability in 
blazars have so far focused on HBLs and, in particular, on the attempt to identify clear pat- 
terns of time lags between hard and soft X-rays. However, such studies have yielded rather 
inconclusive and often contradictory results (e.g., for Mrk 421: Takahashi et al. (1996); Fos- 
sati et al. (2000); Takahashi et al. (2000); or PKS 2155-304: Chiapetti et al. (1999); Zheng 
et al. (1999); Kataoka et al. (2000); Edelson et al. (2001)). Instead, the so-called “spectral 
hysteresis” of blazar X-ray spectral variability may prove to be a more promising diagnostic 
of the physical nature of acceleration and cooling processes in blazar jets: When plotting 
the X-ray spectral hardness vs. the X-ray flux (hardness-intensity diagrams = HIDs), some 
HBLs (e.g., Mrk 421 and PKS 2155-304) have been observed to trace out characteristic, 
clockwise loops (Takahashi et al. (1996); Kataoka et al. (2000)). In terms of pure SSC jet 
models, such spectral hysteresis can be understood as the synchrotron radiation signature of 
gradual injection and/or acceleration of ultrarelativistic electrons into the emitting region, 
and subsequent radiative cooling (Kirk, Rieger, & Mastichiadis 1998; Georganopoulos Sc 
Marscher 1998; Kataoka et al. 2000; Kusunose, Takahara, & Li 2000; Li Sc Kusunose 2000). 
However, interestingly, such spectral hysteresis could not be confirmed in a recent series of 
XMM-Newton observations of Mrk 421 (Sembay et al. 2002). 

In LBLs, the soft X-ray emission is also sometimes dominated by the high-energy end of 
the synchrotron emission component (Tagliaferri et al. 2000; Ravasio et al. 2002), so similar 
spectral hysteresis phenomena should in principle be observable. However, those objects 
are generally much fainter at X-ray energies than their high-frequency peaked counterparts, 
making the extraction of time-dependent spectral information an observationally very chal- 
lenging task (see, e.g., Bottcher et al. (2002)), which may require the new generation of X-ray 
telescopes such as Chandra or XMM-Newton. Extracting the physical information contained 
in the rich X-ray variability patterns exhibited by BL Lac objects requires detailed theoret- 



ical modeling of the time-dependent particle acceleration and radiation transport processes 
in the jets of blazars. Previous analyses of these processes (Kirk, Rieger, & Mastichiadis 
1998; Georganopoulos & Marscher 1998; Chiaberge &c Ghisellini 1999; Kataoka et al. 2000; 
Kusunose, Takahara, & Li 2000; Li & Kusunose 2000; Krawczynski, Coppi, & Aharonian 
2002) have led to significant progress in our understanding of the particle acceleration and 
radiation mechanisms in HBLs, but were restricted to pure SSC models, with parameter 
choices specifically targeted towards HBLs, so those results may not be directly applicable 
to intermediate or low-frequency peaked BL Lac objects or even FSRQs. A notable exception 
is a recent study by Sikora et al. (2001), who included a significant contribution of external 
Compton radiation to the high-energy emission of blazars, and focused on the modeling of 
photon-energy dependent light curves and time lags between different frequency bands. They 
applied their results to the FSRQ 3C 279, and concluded that the correlated X-ray/ 7-ray 
variability of this quasar was inconsistent with X-rays and 7-rays being produced by the 
same radiation mechanism because otherwise significant systematic time lags between the 
7-ray and X-ray flaring behaviour would be expected, contrary to the observations (e.g., 
Hartman et al. (2001)). 

In the present paper, we describe a newly developed combined SSC + ERC jet radiation 
transfer code, accounting for time-dependent particle acceleration and injection, radiative 
cooling, and escape, coupled to the self-consistent treatment of the relevant photon emission, 
absorption, and escape processes. In §2 we give a brief description of the underlying blazar jet 
model. The numerical procedure used in our code will be outlined in §3. We present results 
of a detailed parameter study, relevant for application to intermediate and low-frequency 
peaked BL Lac objects, in §4. We summarize in §5. 


2. Model Description 

The blazar model used for this study is a generic leptonic jet model. It is assumed 
that a population of ultrarelativistic, non-thermal electrons (and positrons) is injected at a 
generally time-dependent rate into a spherical emitting volume of co-moving radius R b (the 
“blob”). The injected pair population is specified through an injection power L inj (t) and 
the spectral characteristics of the injected non-thermal electron distribution. We assume 
that electrons are injected with a single power-law distribution with low and high energy 
cutoffs 71 and 72, respectively, and a spectral index q so that the injection function Q e ( 7; t ) 
[cm -3 s -1 ], in the co-moving frame of the emitting region, is 


with 
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where is the blob volume in the co-moving frame. 
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The jet is powered by accretion of material onto a supermassive central object, which 
is accompanied by the formation of an accretion disk which, for the purpose of this study, 
we have represented by a standard Shakura-Sunyaev disk with a bolometric luminosity of 
Ld = 10 45 ergs s -1 . The choice of this and several other standard parameters is motivated by 
a recent modeling study of the LBL W Comae (Bottcher, Mukherjee, k Reimer 2002). The 
randomly oriented magnetic field B is fixed through the choice of an equipartition parameter 
e B , which determines the fraction of the magnetic field energy density u B compared to its 
value for equipartition with the relativistic electron population in the emission region. The 
blob moves with relativistic speed v/c = fi T = ^/l — 1/P along the jet which is directed 
at an angle 0 o bs (with // = cos0 O b s ) with respect to the line of sight. The Doppler boosting 
of emission from the co-moving to the observer’s frame is determined by the Doppler factor 


As the emission region moves outward along the jet, particles are continuously being 
injected, are cooling, primarily due to radiative losses, and may leak out of the system. We 
parametrize particle escape through an energy-independent escape time scale t e sc = r] R^/c 
with T) > 1. Radiation mechanisms included in our simulations are synchrotron emission, 
Compton upscattering of synchrotron photons (SSC = Synchrotron Self Compton scat- 
tering), and Compton upscattering of external photons (EC = External Compton scat- 
tering), including photons coming directly from the disk as well as re-processed photons 
from the broad line region. The broad line region is modelled as a spherical shell between 
r BLR,in = 0.2 pc and t blr, out = 0.25 pc, and a radial Thomson depth tt,blr which is consid- 
ered a free parameter. 77 absorption and the corresponding pair production rates are taken 
into account self-consistently. Motivated by the ~ 10 hr minimum variability time scale 
observed in W Comae (Tagliaferri et al. 2000), we choose Rb ~ 10 16 cm, and T = D = 10, 
which implies 6 ohs = 5.74°. 


Based on an equipartition parameter e B ~ 1, we expect typical magnetic field values 
of order B ~ 1 G (Bottcher, Mukherjee, k Reimer 2002), which implies a synchrotron 
cooling time scale (in the observer’s frame) of electrons emitting synchrotron radiation at an 
observed energy E sy = 1 £' keV keV of 


E^ 2 hr. 
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which, for X-ray photon energies, is shorter than the dynamical time scale R B /(Dc), in 
agreement with the approximately symmetric shape of the X-ray light curves generally ob- 
served both in W Comae (Tagliaferri et al. 2000) and BL Lacertae (Ravasio et al. 2002). 


3. Numerical Procedure 

In order to treat the time-dependent electron dynamics and radiation transfer problem 
in the emitting volume, we solve simultaneously the kinetic equation for the relativistic 
electrons, 


dnei 7,0 

dt 

and for the photons, 
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Here, {d , y/dt)\ oss is the radiative energy loss rate for the electrons, Qe{7,t) is the sum of the 
external injection rate Q ^ from Eq. 1 and the intrinsic 77 pair production rate, n p h,em(c0 
and n p habs(c^) are the photon emission and absorption rates corresponding to the various 
radiation mechanisms, and * ph , esc = (3/4) R b /c. In Eq. 4, electron cooling is approximated 
as a continuous function of time (i.e. , the energy of an individual electron is described as a 
differentiable function of time). This would be inaccurate if a significant contribution to the 
cooling rate were due to Compton scattering in the Klein-Nishina limit since in that case, the 
electron is transferring virtually all of its energy to a soft photon in a single scattering event. 
However, in the parameter ranges which we are primarily interested in, electron cooling is 
dominated by synchrotron losses and Compton scattering in the Thomson regime, for which 
Eq. 4 is a good approximation. In Eq. 5, the emissivity term contains the contribution 
from Compton scattering into a given photon energy interval. Since in all model situations 
considered here, the Thomson depth of the emitting region is 7 t “C 1CT 6 , the modification 
of the photon spectrum due to scattering of photons out of a given photon energy range is 
negligible. 

The relevant electron cooling rates and photon emissivities and opacities are evaluated 
using the well-tested subroutines of the radiation jet transfer code described in detail in 
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Bottcher, Mause, & Schlickeiser (1997) and Bottcher & Bloom (2000). The full Klein-Nishina 
cross section for Compton scattering and the complete, analytical solution for the 77 pair 
production spectrum of Bottcher & Schlickeiser (1997) are used. The discretized electron 
continuity equation can be written in the form of a tri-diagonal matrix as in Chiaberge & 
Ghisellini (1999), which can be readily solved using the standard routine of Press et al. 
(1992). This procedure turns out to be very stable if, instead of the sharp cutoffs of the 
electron injection function (1), we introduce continuous transitions to very steep power-laws 
to mimic these cutoffs. Specifically, we add a low-energy branch with Q e (T,t ) ex j 2 for 
7 < 71, and <2e(7;0 oc 7~^+ 3 ) for 7 > 72. After each electron time step, we update the 
photon distribution using a simple explicit forward integration of the discretized Eq. 5. 

We have carefully tested our code by running it with parameters identical to those used 
for Figs. 6 - 14 of Li & Kusunose (2000), who are using a very similar numerical approach. 
We find very good agreement with their results, with only minor discrepancies which are 
due to our replacing the high-and low-energy cutoff of the electron injection spectrum by 
continuous transitions to very steep power-laws as described above. Specifically, this results 
in more moderate spectral indices of the synchrotron spectra just beyond the high-energy 
cutoffs. 


4. Numerical Results 

We have performed a large number of simulations, studying the influence of various 
model parameters on the resulting broadband spectra, light curves, and X-ray hardness- 
intensity diagram tracks. In each one of our simulations, we have assumed an underlying, 
quiescent injection power of L|' n “ = 10 38 ergs s _1 , on top of which we inject particles with 
various flaring injection powers in the range 10 4 ° ergs s _1 < Z,f nj < 10 43 ergs s _1 . As a 
standard model setup, we choose an injection electron function given by 7 1 = 10 3 , 72 = 10 5 , 
and q = 2.5. In the base model, we have L-L = 10 41 ergs s _1 , extending as a step function 
in time over 2 dynamical time scales, A ffo = 2 R^/c in the co-moving frame. The injection 
event is centered around a distance %o = 0.1 pc from the central engine. The BLR Thomson 
depth is chosen to be 0 in the base model. 

Subsequently, we investigate the influence of changing (1) the flaring injection power, 
(2) the BLR Thomson depth and, accordingly, the contribution of external photons to the 
soft photon field for Compton scattering, (3) the electron injection spectral index q, (4) the 
duration of the flaring injection event, (5) the time profile of the flaring injection event, (6) 
the electron escape time scale parameter. The parameters used for the individual runs are 
quoted in Tab. 1. 
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Our base model is simulation no. 2. In Figs. 1 and 2, we have compiled a sequence 
of co-moving electron spectra, snap-shot SEDs, and the time averaged photon spectrum 
resulting from our base model. Light curves at 3 selected X-ray energies as well as in the 
optical (R-band) and at hard X-rays (30 keV) are plotted in Fig. 3, and tracks in the 
hardness-intensity diagrams (HIDs) at three different X-ray energies are compiled in Fig. 
4. The figures illustrate the gradual build-up of the electron density in the emission region, 
competing with radiative cooling, which is faster than the injection time scale at electron 
energies of 7 ^ 10 4 . Radiative cooling is the dominant process affecting the electron spectra 
after the end of the flaring injection episode at t = 2 f dyn . The time-dependent photon spectra 
as well as the light curves demonstrate that we do not expect significant peak time delays 
within the synchrotron component at frequencies u > 10 14 Hz, but that the high-energy 
(SSC) component is delayed by 1 dynamical time scale due to the gradual accumulation 
of seed photons for Compton scattering. The figure also indicates the very moderate flux 
variability at energies just above the synchrotron cut-off, which is located at ~ 1 keV in 
our example. Fig. 4 illustrates the spectral hysteresis phenomenon (keeping in mind that 
additional contributions from previous injection episodes should close the tracks in the sense 
that they are expected to start out near the end points of the tracks shown in the figure). 
In agreement with Li & Kusunose (2000) we find that — at least in this generic case the 
spectral hysteresis tracks can change their orientation from clockwise to counterclockwise 
as one goes from photon energies below the synchrotron cutoff to energies above the cutoff, 
where the spectrum is dominated by Compton scattering (SSC). 

In the following, we are focusing on the time-averaged photon spectra, the light curves, 
and the X-ray spectral hysteresis, and investigate how those aspects are affected by variations 
of individual parameters. 


4.1. Electron Injection Power 

The effect of an increasing injection power — corresponding to a higher density of 
injected, relativistic particles in the emitting region — is illustrated in Figs. 5-7. In 
addition to a corresponding increase in the overall bolometric luminosity, this leads also 
to a stronger relative energy output in the SSC-dominated Compton emission at X- and 
7-ray energies, as expressed, e.g., in Eq. (19) of Chiang Bottcher (2002). The photon 
spectral index of the time-averaged emission at optical - soft X-ray frequencies remains 
robust at a 0 _ x ~ 1.25 due to optically thin synchrotron emission from the cooled electron 
spectrum with injection spectral index q = 2.5. As the bolometric luminosity (and the 
electron cooling) becomes dominated by the SSC mechanism, one would expect that this 
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changes to the canonical value of a 0 _ x = 1.5, which is a result of the decaying electron 
cooling rate in an SSC-dominated cooling scenario (Chiang k Bottcher 2002). However, in 
the situation simulated here, the first-order SSC peak is rapidly (within ~ 3t dyn ) decaying 
into the keV energy range and dominating over the instantaneous synchrotron emission at 
UV - X-ray energies. This leads to a significant hardening of the time-averaged optical - 
X-ray spectrum, which even becomes inverted in i/F v space in our most extreme test case 
(simulation no. 8). 

In the light curves (see Fig. 6), the more rapid electron cooling with increasing electron 
injection power manifests itself in an overall increasing amplitude of variability at all energies. 
In particular, as SSC cooling becomes more important, even the harder X-rays begin to 
exhibit significant variability on the dynamical time scale, in contrast to the synchrotron- 
cooling dominated cases. Furthermore, while for very low injection powers, the synchrotron 
cooling time scale for optical synchrotron emission is comparable to the injection time scale, 
resulting in a time delay of ~ a few hr between X-ray and optical emission, the optical light 
curve peaks at the end of the injection episode for higher injection powers, simultaneously 
with the X-rays. However, when SSC cooling becomes dominant, the gradually increasing 
energy density in the soft photon field during the injection episode actually has the effect 
that the X-ray light curves are peaking at the beginning of the injection episode, which 
would, again, lead to a time delay of ~ a few hr between X-ray and optical flares. 

Fig. 7 illustrates how the tracks in the HIDs at different X-ray energies are drastically 
changing for different injection powers. In particular at X-ray energies just below or at the 
synchrotron cutoff (~ 1 keV), the flux maxima are occurring at significantly different values 
of the local spectral index a for different values of the injection power. Specifically, the local 
spectral indices at the time of the peak peak flux are significantly smaller (harder) for larger 
values of the injection power. Obvious changes in the orientation of the spectral hysteresis 
tracks are not found in these simulations. 

As mentioned earlier (see Eq. 3), in the test cases investigated here, the synchrotron 
cooling time scale of electrons emitting synchrotron radiation at X-ray energies, is shorter 
than the dynamical time scale, which is of the same order as the injection time scale. Conse- 
quently, our results can be qualitatively compared to those of Li & Kusunose (2000) for the 
“short cooling time limit”, bearing in mind that the parameter values in our simulations have 
been chosen appropriate for intermediate and LBLs, while Li k Kusunose (2000) focused on 
the application to the HBL Mrk 421. In particular the light curves displayed in their Figs. 
9-11 exhibit the same general trends as we have found in our set of simulations. 
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4.2. External Photons 

In order to investigate the influence of an increasing contribution of external photons 
to the soft photon field for Compton scattering, we performed a series of simulations with 
increasing values of t TjBL r, fr° m 0 to 1. The resulting time-averaged photon spectra are 
shown in Fig. 8, which clearly shows the emergence of the external Compton (EC) com- 
ponent at GeV 7-ray energies. The impact of this additional emission component on the 
lower-frequency emission is small as long as its bolometric energy output is smaller than or 
comparable to the synchrotron vF„ flux. Only when EC cooling becomes dominant over 
synchrotron cooling, are the effects on the synchrotron + SSC dominated portion of the 
spectrum (radio frequencies - MeV 7-rays) noticeable. This is the case when the comoving 
energy densities of the magnetic field and the external photons, Ug and u ext , respectively, 
become comparable. From 


u, 


ext 


Ut 


L45 TT,BLR r? 

(B/G) 


2o "/r,i’rrv - (6) 

— where L 45 is the accretion disk luminosity in units of 10 45 ergs s -1 , IT = T/IO, and 
r 0 . 2 is the radius of the inner boundary of the BLR in units of 0.2 pc — we can estimate 
that this happens at 7t,blr ~ 0.1. Specifically, for tt,blr > 0.1, the effect of dominant 
external-Compton cooling results in a reduction of the time-averaged emission around the 
synchrotron peak — leading to a spectral hardening of the synchrotron emission — , and a 
suppression of the SSC emission. The suppression of the low-frequency synchrotron emission 
can be explained as the combined effect of two causes: First, electrons at energies above the 
low-energy cut-off are radiatively cooling on a time scale much shorter than the dynamical 
one (see Eq. 7 below). Consequently, the particle spectrum of high-energy electrons i njected 
during the flare is rapidly depleted within less than one dynamical time scale from the end 
of the flaring episode. However, after this episode, we are still injecting electrons (although 
at a much smaller rate corresponding to L?“) into the blob, which continue to present a 
high-energy electron population. In the case of the extremely fast cooling rate at t t ,blr ~ T 
this additional high-energy electron population leads to a significant additional contribution 
of synchrotron emission at intermediate energies beyond the down-shifted maximum energy 
of electrons injected during the flare and, consequently, to a hardening of the synchrotron 
spectrum. 


The second cause of spectral hardening of the low-frequency synchrotron spectrum is 
related to our parametrization of the magnetic field in terms of an equipartition parameter 
For electron injection spectral indices q > 2, most of the energy in the electron population 
is carried by electrons near the low-energy cutoff. Consequently, the co-moving magnetic 
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field will decay on a time scale given by the radiative cooling time scale of the lowest-energy 
electrons, which is 


4 = ^ec( 7 i) ~ 4 x 10 4 


0.2 


r T,BLR 7l,3 L 45 Tf 


(7) 


where 7i,3 = 7i/10 3 . For r T)BLR > 0.1, the time scale (7) is of the order of or shorter than the 
dvnamical time scale. In contrast, the energy density of the external radiation field remains 
approximately constant during the evolution of the blob. Consequently, as radiative cooling 
proceeds and shifts the electron distribution towards lower energies, a steadily decreasing 
fraction of the electron energy will be converted to synchrotron radiation. This leads to a 
spectral hardening of the time-averaged synchrotron spectrum compared to the canonical 
a — q/2 spectrum above the break frequency 17 = i/ L 0 of a cooling electron population 
in a constant magnetic field B = B 0 with u Lfl = eB 0 /(2nm e c). At frequencies below u u 
the time-averaged spectrum would have a slope of a = 1/2 in the constant-magnetic-field 
case. The steepening of this low-frequency part of the synchrotron spectrum in the case of 

a magnetic field proportional to the equipartition value can be derived analytically in the 
following way. 


At late times, flaring electron distribution has basically collapsed to a 5 function in 

electron energy, i. e. n(j e ,t) oc S(y e - 'y[t]). Then, the synchrotron emission coefficient will 
behave as 


3v,sy(t) OC % y S(w - l[t) 2 V L [t)) (8) 

oc l[t] z 5(v -^[tf /2 v Lfi ). (9) 

where j ay oc 7 (t) 3 is the synchrotron loss rate. The additional factor of 7 results from the 
equipartition prescription, B 2 /8n = The electron evolution is governed by EC losses, 

and since the EC photon field is essentially constant over the flare episode, the electron 
Lorentz factor still evolves according to 7 » -(4/3)(cr T /m e c)u' xt j 2 . The time-integrated 
synchrotron spectrum is then 


oc 

where we have used z. oc 7 5 / 2 . The short-dashed curve in Fig. 8 shows the result of a test 
calculation in which we held the magnetic field constant, while all other parameters were 


J dtj[t} 3 5(v - 7[£] 5/ V L)0 ) 

r ,3 -1 

7W — oc v 


dt 


-1/5 


( 10 ) 

( 11 ) 
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identical to the t T)B lr = 1 simulation, in order to verify that the spectral hardening at radio 
frequencies is indeed partially a consequence of our magnetic field parametrization. 

Fig. 9 shows that the impact of a strong external Compton component on the optical 
and X-ray light curves is very moderate. In particular, the light curves at energies below 
the synchrotron peak remain virtually unchanged. However, a strong external Compton 
component, dominating the bolometric luminosity (t t ,blr ^ 0.1 in our case) leads to a 
significantly faster decay of the light curves at X-ray energies beyond the synchrotron peak. 

Just as the light curves, also the X-ray spectral hysteresis characteristics remain virtually 
unchanged, even in the case of a strongly dominant external Compton component, except 
for a moderate softening of the local spectrum in the decaying phase of the flare at energies 
beyond the synchrotron peak (see Fig. 10). 


4.3. Electron Spectral Index 

The value of the electron injection spectral index q should be rather easily determined by 
measuring the time-averaged optical - X-ray spectral index of the strongly cooled synchrotron 
spectrum. As illustrated in Fig. 11, this spectral change is accompanied by a shift of the 
SSC peak towards higher frequencies as the injection spectrum hardens: For q > 2.5, the 
SSC peak is located at e S sc ~ 7?e sy , while for q < 2.5, it shifts towards e S sc ~ 7 As 
illustrated in Fig. 12, the characteristics of the light curves are only marginally different for 
different values of q. The X-ray spectral hysteresis at X-ray energies below the synchrotron 
cutoff is obviously shifted according to the change in injection spectral index, but its basic 
characteristics remain unchanged (see Fig. 13). An interesting qualitative change of the 
spectral hysteresis can be seen for energies just above the synchrotron cutoff: While for hard 
electron injection spectra (q < 2.5), the peak flux is reached at a steep local X-ray spectrum 
(i.e. the spectrum is dominated by the synchrotron component at the time of peak flux), a 
soft injection spectrum (q > 2.5) leads to a hard local spectral index at the time of peak flux 
(i.e. the spectrum is dominated by the SSC component at that time). 


4.4. Duration of the Flare 

In order to investigate the influence of the duration of the electron injection event 
causing the flare, we have performed simulations with 3 different values of t m j, where we 
kept the total energy input during the flare constant (simulations no. 2, 11, and 12). The 
time-averaged SEDs from those simulations are virtually indistinguishable. Fig. 14 shows 



that the longer injection time scale at a lower injection power leads to a more gradual rise 
of the light curves at energies above the synchrotron cut-off, whereas at lower frequencies, 
this leads to a rather marginal modification of the rising portion of the light curve, followed 
by an extended plateau until the end of the injection episode. The decaying portion of 
the light curves at optical and X-ray frequencies seems to be virtually independent of the 
duration of the injection event. In the X-ray spectral hysteresis (Fig. 15), we find that 
for longer electron injection events, the rising-flux portion of the track in the HID at X-ray 
energies below the synchrotron cutoff occurs with increasingly softer local spectra, whereas 
the decaying portion remains virtually unchanged. At energies above the synchrotron peak, 
the trend concerning the rising portion is opposite: With longer duration of the injection 
episode, the local spectra are becoming harder. 

From an observational point of view, one should be able to distinguish situations in 
which the duration of the flaring event is substantially longer than the dynamical time scale, 
by virtue of the extended high-flux plateaus at energies below the synchrotron peak, which 
are generally not observed in BL Lac objects in which there is evidence for a substantial 
contribution from synchrotron emission to the X-ray flux (e.g., Tagliaferri et al. (2000); 
Ravasio et al. (2002)). This seems to indicate that the generic situation of t in j ~ 2 i dyn might 
be a realistic assumption. 


4.5. Time Profile of the Flare 

The step function time profile of the electron injection power during the flare is certainly 
a rather crude over-simplification of any realistic acceleration scenario. In order to investigate 
whether this particular choice of the time profile has a significant impact on our results, we 
have calculated an additional set of simulations (nos. 13 - 15), with triangular injection 
profiles. Here, we have introduced a linear rise and decay of the injection power on time 
scales t r and tj, respectively, and have chosen the maximum of the profile to be twice the 
injection power of the step-function case in order to keep the total injected energy at the 
same value as in our base model. 

The time- averaged photon spectra for all of these cases are virtually identical. Equally, 
the light curves at energies above the synchrotron cutoff are only marginally affected by the 
details of the injection time profile (see Fig. 16), while at lower energies the light curves tend 
to track the injection time profile to a certain extent during the rising portion of the light 
curves. The decaying portion of the light curves is generally independent of the injection time 
profile for all optical and X-ray photon energies. Fig. 17 illustrates that the impact of the 
detailed injection profile on the X-ray spectral hysteresis characteristics is rather moderate. 
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even for the extreme (and very artificial) cases of the time profiles c (gradual rise over 2 f dyn , 
sharp decay) and d (sharp rise, gradual decay over 2 f d yn ) illustrated here. Its influence is 
basically restricted to the flux-rise portion of the HID track and to photon energies below the 
synchrotron cut-off, where the local spectra tend to be harder for time profiles with maxima 
closer to the onset of the flare. The HID tracks for all our test cases show almost identical 
spectral indices at the time of maximum flux. 


4.6. Electron Escape Time Scale 

At the low-energy end of the electron spectra, particles cooling down from higher energies 
will either accumulate and build up a 7 2 power-law spectrum at energies below 71 > escape, 

depending on the value of the escape time scale parameter 77 . We can define a critical escape 
parameter 


ricr = 


= 2.3 


(TV 

\ 10 16 cm ) \G / 


(—) 

VlO 3 / 


-1 


(12) 


RT S y{ 7l) ~~V 1016cn n 

for which the escape time scale for electrons at energy 71 equals the synchrotron cooling time 
scale. For 77 < r] CT , the electron distribution will maintain a sharp low-energy cutoff at 71 , 
while for 77 r] CT , a -2 low-energy power law will develop. In order to investigate whether 
our choice of 77 = 10 has a significant impact on our results, we have done test simulations 
with 77 = 3 and 77 = 30, respectively. We find all relevant aspects — the time-averaged 
spectra, the monoenergetic light curves, and the spectral hysteresis curves to be virtually 
independent of 77 within reasonable bounds. 


4.7. Other parameters 

In the previous subsections, we have discussed the impact of various parameters on the 
broadband SED and the X-ray variability characteristics of models for intermediate and low- 
frequency peaked BL Lac objects. There are still a few more parameters left which we have 
fixed in our suite of test simulations. In particular, the choice of the cutoffs of the electron 
distribution, 71 and 72 , the magnetic-field equipartition parameter eg, the Doppler factor 
D, and the size of the emitting region, Rb may have an impact on the spectral variability 
characteristics. However, those parameters can generally be constrained rather well through 
the observed overall spectral characteristics, and through variability time scale considerations 
(see, e.g., Tavecchio, Maraschi, & Ghisellini (1998); Costamante &; Ghisellini (2002)). The 
values adopted here are representative for typical LBLs like BL Lacertae (Madejski et al. 
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1999 ; Bottcher & Bloom 2000) or W Comae (Tagliaferri et al. 2000; Bdttcher, Mukherjee 
& Reimer 2002). Furthermore, significantly different values of 72, and D would shift the 
synchrotron cutoff out of the X-ray regime so that the diagnostics developed here may not 
be applicable to the X-ray variability of BL Lac objects. 


5. Summary and Conclusions 

We have presented a detailed parameter study of the time-dependent electron injection 
and kinematics and the self-consistent radiation transport in jets of intermediate and low- 
frequency peaked BL Lac objects. Those objects are currently of great interest as the 
steadily improving capabilities of current and future air Cerenkov detector facilities might 
allow the detection of this class of blazars at multi-GeV energies in the near future. At 
the same time, some of these objects exhibit interesting X-ray variability features which can 
now be studied in detail with the new generation of X-ray telescopes, in particular Chandra 
and XMM-Newton. Furthermore, the GLAST mission, scheduled for launch in 2006 , is 
expected to detect many more BL Lac objects at multi-MeV - GeV energies and bridge the 
observational gap between the energy ranges previously covered by the EGRET instrument 
on board CGRO, and the ground-based air Cerenkov facilities. 

In our study, we have focused on the impact of various specific parameter choices and 
variations on the broadband SEDs, optical and X-ray light curves, and the spectral hysteresis 
phenomena previously observed in several high-frequency peaked BL Lac objects, but also 
expected to be observable in intermediate and low-frequency peaked BL Lacs. 

Very important conclusions can be drawn from a comparison of our results concerning 7- 
ray bright sources dominated by either SSC or external Compton emission. For a given level 
of flux at GeV energies, at a level comprable to or exceeding the synchrotron vF v peak flux, 
those two scenarios should be clearly distinguishable by virtue of the X-ray variability during 
flaring episodes: If roughly symmetric flare time profiles at soft X-rays below the synchrotron 
peak are observed, and the local, time-resolved X-ray spectra are soft - consistent with a 
case with negligible 7-ray emission — the 7-ray emission might be dominated by external- 
Compton emission. However, if the X-ray time profiles show a clear sign of a very rapid 
rise and more gradual decay, and the time-resolved spectra are significantly harder than 
corresponding to a case without strong 7-ray emission, we expect a strong contribution from 
the SSC mechanism to the 7-ray emission. Most notably, this diagnostic of the 7-ray emission 
mechanism does not require a detailed spectral measurement — which will be hard to achieve 
and require long integration times, even with GLAST -, but only a rough estimate of the 
GeV flux. It relies primarily on X-ray variability studies, combined with the type of detailed 
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time-dependent radiation modeling presented in this paper. 

We have shown that our results are not significantly impacted by the special choice 
of poorly determined parameters like the details (exact duration and time profile) of the 
acceleration / injection events leading to flaring activity, or the electron escape time scale 
parameter. Consequently, the X-ray variability of the high-frequency end of the synchrotron 
emission in intermediate and low-frequency peaked BL Lac objects can be used as a very 
robust diagnostic to unveil the nature of the high-energy emission in this type of blazars. 

The work of MB is supported by NASA through Chandra Postdoctoral Fellowship Award 
no. 9-10007, issued by the Chandra X-ray Center, which is operated by the Smithsonian 
Astrophysical Observatory for and on behalf of NASA under contract NAS 8-39073. 
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Fig- 1. — Sequence of co-moving electron spectra in the emission region for our base model, 
simulation no. 2. The curves are labeled by time in multiples of the dynamical time scale, 
^dyn = Rb/(Dc) = 3.33 x 10 4 s. For model parameters, see Tab. 1. 
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Fig. 2. — Sequence of spectral energy distributions (upper panel) and the time averaged SED 
(lower panel) from our base model, simulation no. 2. The dotted vertical lines indicate the 
frequencies at which light curves and hardness-intensity correlations have been extracted; 
the dashed vertical lines indicate the remaining two frequencies at which light curves have 
been extracted. 
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Fig. 3. Light curves at three different X-ray energies, in the optical (R-band) and at 
hard X-rays for our base model. Upper panel: Absolute flux values; lower panel: fluxes 
normalized to their peak values. The vertical dotted line indicates the end of the flaring 
injection episode. 
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Fig. 4. — Tracks of the simulated spectra in the hardness-intensity diagrams at three different 
X-ray energies for our base model, a is the local energy spectral index F„ a v~ a at the 
respective photon energy. Stars indicate the locations at multiples of the dynamical time 
scale during the simulation. 
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Fig. 5.-— Time averaged photon spectra for different values of the flaring injection power 
from simulations no. 1, 2, 3, and 8. The dotted vertical lines indicate the frequencies at which 
light curves and hardness-intensity correlations have been extracted; the dashed vertical lines 
indicate the remaining two frequencies at which light curves have been extracted. 


i TTIl 1 — i i > mil 




10 

10 

10 

V 

10 



w no 5 «] 


30 keV 


6. — Optical and X-ray light curves for different values of the flaring injection power 
from simulations no. 1, 2, and 3. 







Fig- 7. — Tracks in the harness-intensity diagrams at X-ray energies for different values of 
the flaring injection power Lf nj , from simulations no. 1, 2, 3, and 8. Symbols are the same as 
in Fig. 5: short-dashed (10 40 erg/s), solid (10 41 erg/s), dot-dashed (10 42 erg/s), long-dashed 
(10 43 erg/s). Stars indicate the locations at multiples of the dynamical time scale during the 
respective simulations. 
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Fig. 8. — Time averaged photon spectra for different intensities of the external soft radiation 
field, parametrized through different values of 7 t,blr. 7 from simulations no. 2, 5, 6, 7, and 18. 
The short-dashed curve shows a test simulation with constant magnetic field, while all other 
parameters were identical to the tt,blr = 1 simulation, illustrating the effect of our magnetic 
field parametrization as a constant fraction of the equipartition magnetic field. The dotted 
vertical lines indicate the frequencies at which light curves and hardness-intensity correlations 
have been extracted; the dashed vertical lines indicate the remaining two frequencies at which 
light curves have been extracted. 
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Fig. 9. — Optical and X-ray light curves for different intensities of the external soft radiation 
field, parametrized through different values of t^blr, from simulations no. 2, 6, and 7. 
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Fig. 10 — Tracks in the harness-intensity diagrams at X-ray energies for different intensities 
of the external soft radiation field, parametrized through different values of tt.blr, from 
simulations no. 2, 6, and 7: solid (t T)B lr = 0), dashed (tt,blr = 0.1), dot-dashed (tt,blr = 
1). Stars indicate the locations at multiples of the dynamical time scale during the respective 
simulations. 
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Fig. 11. — Time averaged photon spectra for different values of the electron injection spec- 
tral index q , from simulations no. 2, 9, and 10. The dotted vertical lines indicate the 
frequencies at which light curves and hardness-intensity correlations have been extracted; 
the dashed vertical lines indicate the remaining two frequencies at which light curves have 
been extracted. 
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Fig. 12. — Optical and X-ray light curves for different values of the electron injection spectral 
indes q, from simulations no. 2, 9, and 10. 
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Fig. 13. — Tracks in the harness-intensity diagrams at X-ray energies for different values of 
the electron injection spectral index q , from simulations no. 2, 9, and 10: dashed ( q = 2.2), 
solid ( q = 2.5), dot-dashed (q = 2.8). Stars indicate the locations at multiples of the 
dynamical time scale during the respective simulations. 
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Fig. 14. — Optical and X-ray light curves for different values of the duration of the flar- 
ing injection time scale, keeping the total energy input during the injection event constant 
(simulations no. 2, 11, and 12). 
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Fig. 15. — Tracks in the harness-intensity diagrams at X-ray energies for different values 
of the flaring injection time scale, keeping the total energy input during the injection event 
constant (simulations no. 2, 11, and 12): solid (t inj = 2 t dyn ), dashed (t inj = 4f dyn ), dot- 
dashed (<j n j = 6 f d yn ) • Stars indicate the locations at multiples of the dynamical time scale 
during the respective simulations. 
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Fig. 16. — Optical and X-ray light curves for different time profiles of the electron injection 
power, keeping the total energy input during the injection event constant (simulations no. 
2, 14, and 15). The small insets illustrate the injection time profiles (L in] vs. time). 
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Fig. 17. — Tracks in the harness-intensity diagrams at X-ray energies for different time 
profiles of the electron injection power, keeping the total energy input during the injection 
event constant (simulations no. 2, 14, and 15): solid (profile a = step function), dashed 
(profile c = triangular profile with t r = 2f dyn and t d = 0), dot-dashed (profile d = triangular 
profile with t r = 0 and t d = 2 f d yn ). Stars indicate the locations at multiples of the dynamical 
time scale during the respective simulations. 
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Table 1 . Parameters of the simulations used for our parameter study. Throughout the 
series of simulations, we have used L?" = 10 38 ergs s _1 (electron injection luminosity during 
quiescence), 71 = 10 3 (low-energy cutoff of injected electron spectrum), 72 = 10 5 
(high-energy cutoff of injected electron spectrum), D = 10 (Doppler boosting factor), 

R b - 10 16 cm (blob radius), and e B - 1 (magnetic- field equipartition parameter). The time 
profiles in the 6th column are: (a) step function, (b) triangular with linear rise and decay 
with equal time scales (t r> d = tinj/2), (c) linear rise and instantaneous drop ( t T = i i n j ) > (d) 
instantaneous rise and linear decay ( t ^ = fj n j)- For time profiles (b) - (d), the parameter 
Lf n j is the maximum injection power, tj is the electron escape time parameter, defined by 
t e ,esc = V Rb/c. The parameters changed with respect to the base model (no. 2 ) printed in 

boldface. 
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*In this simulation, the magnetic field was held constant at B = 0.4 G. 



